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ABSTRACT 


This  is  the  third  interim  report  of  a  study  whose  purpose  is  to 
develop  and  apply  a  generalized  theory  of  coupled  modes  of  propagation 
to  the  study  of  interactions  in  distributed  microwave  devices. 

In  Sec.  I  we  develop  methods  for  the  application  of  group  theoretic 
techniques  to  conservative  linear  systems  that  may  or  may  not  be  uniform. 
It  is  shown  that  transmission  matrices  that  solve  a  particular  problem 
do,  in  fact,  form  a  group.  This  group  is  a  continuous  (Lie)  group  and 
can  be  characterized  by  its  infinitesimal  transformations.  It  is  shown 
that  the  appropriate  infinitesimal  transformations  can  be  obtained  from 
the  system  operator,  Rlz).  We  can  then  investigate  whether  or  not  a 
given  system  is  reducible  in  the  group  theoretic  sense.  If  so,  the 
given  problem  can  be  replaced  by  a  simpler  one.  In  any  case,  the 
properties  of  the  system  can  be  discussed  usefully  in  terms  of  the 
commutation  relations  of  tlie  infinitesimal  transformations. 

In  Sec.  II  we  continue  the  work  begun  in  the  previous  quarter  on 
non-uniform  systems  such  that  the  derivative  of  B  is  expressible  as  a 
commutator  of  R  with  another  matrix  A.  We  showed,  last  quarter,  that 
such  a  system  is  explicitly  soluble  if  A  can  be  chosen  to  be  constant. 

We  show,  here,  that  this  is  a  special  case  of  a  more  general  class  of 
explicitly  souble  non-uniform  systems. 

We  have  been  able  to  show  that  the  matrix  A  describes  the  behavior 
of  the  eigenvectors  of  R.  An  explicit  form  for  A  is  obtained  and  some 
of  its  properties  are  developed. 

There  is  reason  to  hope  that  this  type  of  analysis  will  prove  of 
great  generality  and  power  in  the  study  of  non-uniform  systems.  However, 
a  considerable  amount  of  work  remains  to  be  done  to  establish  the 
method  in  the  generality  that  appears  to  be  possible. 


•  • 
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I  INTRODUCTION 


Tlie  following  work  is  based  on  tlic  tJiouglit  tliat  a  conservation  law 
defines  a  group.  That  is,  if  we  consider  all  |)ossible  systems  that 
exhibit  a  given  conservation  law,  we  find  that  their  transmission  matrices 
form  a  continuous  (Lie)  group  of  the  requisite  dimensionality. 

If  we  are  given  a  vector  differential  equation  that  describes  a 
system  with  a  known  conservation  law,  then  this  equation  and  its  boundary 
condition  determines  a  trajectory  through  the  space  of  group  elements. 
Alternatively,  if  the  group  elements  are  suitably  parametrized,  then  the 
vector  differential  equation  and  its  boundary  conditions  determine  a 
trajectory  through  the  space  of  the  narameters. 

The  transmission  matrices  with  different  boundary  conditions  forma 
continuous  set  of  matrices  that  may  be  embedded  in  a  proper  subgroup  of 
the  whole  group  generated  by  the  conservation  law.  If  so,  then  the 
behavior  of  the  system  can  be  described  within  this  subgroup.  Certain 
of  its  properties  then  can  be  obtained  by  consideration  of  the  subgroup. 

In  addition,  it  may  happen  that  this  subgroup  can  be  considered  as  a 
representation  of  a  group  that  is  simpler  than  the  original  group.  Thus, 
for  example,  we  may  be  aljle  to  show  that  a  given  n-dimensiona  1  system  can 
be  expressed  as  a  representation  of  an  appropriate  two-dimensional  group. 
The  n-dimensional  problem,  then,  can  be  replaced  by  an  exactly  analogous 
two-dimensional  one. 

To  establish  this  mode  of  analysis,  we  will  review  in  some  detail 
the  formulation  of  the  problems  that  concern  us  and  the  necessary  parts 
of  the  theory  of  Lie  groups. 

A.  TIIK  SY.STKM  AND  ITS  CONS KIIVAT ION  I.AW 

We  are  concerned  with  linearized  sysi.ems  that  arc  described  by  a 
vector-differential  er|uatiou  of  the  form 


The  n-dimensiona I  column  vector  x  lists  the  significant  system  variables 
which  are  functions  of  the  single  independent  variable,  z.  for  example, 
in  a  TWT,  X  may  list  the  complex  amplitudes  of  tlie  circuit  voltage  and 
current,  and  the  complex  amplitudes  of  the  RF  component  of  space  charge 
or  longitudinal  beam  current  and  velocity.  If  parametric  interaction 
occurs,  we  can  list  these  variables  at  the  various  frequencies  that  are 
coupled . 

The  matrix  R  is  an  n  ^  n  matrix.  If  R  is  independent  of  z,  the 
system  is  said  to  be  uniform.  If  R  is  a  function  of  z,  it  is  non-uniform. 
It  is  the  latter  case  that  is  of  principal  interest  here,  since  it  is 
these  problems  in  particular  that  need  more  subtle  methods  of  analysis. 

In  what  follows  we  shall  not  assume  that  R  is  constant,  except  as 
noted . 

The  solution  to  Eq.  (1)  can  be  written  as 

x(z)  =  M(z)  x(0)  (2) 

where  M(z)  is  an  n  n  matrix  that  is  a  function  of  z  such  that 

M(0)  =  I  (3) 

where  I  is  the  identity  matrix.  The  constant  vector  x(0)  is  the  boundary 
condition  that  is  assumed  to  be  applied  at  z  =  0. 

If  we  substitute  Eq.  (2)  in  Eq.  (I),  and  require  that  it  be  valid 
for  all  x(0),  then  we  find 


<iM 

dz 


-jm 


(4) 


Equation  (4)  with  boundary  condition  (3),  can  be  taken  as  an 
alternative  description  of  the  system. 

We  are  concerned  with  systems  that  exhibit  a  conservation  law — 
conservation  of  net  real  power  flow,  the  linearized  Manley-Rowe  relations, 
Chu's  kinetic  power  theorem,  etc.  .Sucli  a  conservation  law  can  be  expressed 
as  stating  the  invariance  of  a  quadratic  form.  That  is,  we  can  define  a 
constant  matrix  K  such  that  the  scalar,  s,  given  by 


2 


5 


(5) 


=  x*Kx 

is  independent  of  2  if  x  is  a  solution  of  Eq.  (1).  (The  f  means  the 
hermitian  conjugate,  or  complex  conjugate  transpose.  Hence  X^  is  a  row 
vector  and  s  a  scalar.)  We  require  that  K  be  nonsingular,  hermitian 
(K  =  K  +  )  and  independent  of  z. 

It  is  easily  seen,  by  differentiating  Eq.  (5)  and  substituting 
Eq.  (1),  that  the  necessary  and  sufficient  condition  on  R  is  that 

KR  •>  R^K  .  (6) 

We  say  that  R  is  “K-hermitian . ” 

If  we  substitute  Eq.  (2)  in  Eq.  (S)  and  use  the  boundary  condition, 

Eq.  (3),  we  find  that  the  necessary  and  sufficient  condition  on  M  is  that 

M^KM  -  K  (7) 

We  call  such  a  matrix  “/(-unitary.  ” 

In  what  follows,  then,  we  shall  be  concerned  with  the  class  of  all 
/(-unitary  matrices.  We  shall  show  that  this  is  a  group,  which  we  shall 
call  the  “/(-group.” 

As  mentioned  in  the  introduction,  we  consider  the  system  equation, 

Eq.  (4)  and  Eq.  (3)  as  determining  a  trajectory  through  the  designated 
Jf-group. 

B.  THE  /(-GROUP 

We  must,  first,  show  that  Eq.  (7)  does,  in  fact,  determine  a  group. 

A  group  is  defined  in  terms  of  a  set  of  elements  and  a  rule  of  combination. 
In  this  case,  the  set  is  the  set  of  all  matrices,  M,  satisfying  Eq.  (7), 
and  the  operation  is  matrix  multiplication.  For  the  set,  then,  to  form 
a  group,  (1)  it  must  contain  the  identity,  I,  (2)  every  element  must  have 
an  inverse  in  the  set,  and  (3)  if  A  and  B  are  in  the  set,  so  is  AB. 

Clearly  M  ”  I  satisfies  Eq.  (7).  Hence  the  set  contains  the 
identity . 
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Since  K  is  to  be  assumed  nonsingular,  its  determinant  does  not 
vanish.  Since  tlie  determinant  of  the  product  of  matrices  is  the  product 
of  their  determinants,  it  follows  from  Eq.  (7)  that  the  determinant  of 
M  cannot  vanish.  Hence  M” *  exists  and  so  does  If,  now,  we  pre¬ 
multiply  Eq .  (7)  by  and  postmultiply  by  M~*,  we  find  that 

K  =  (8) 

Hence  is  a  member  of  the  set.  Hence,  every  member  of  the  set  has  an 
inverse  that  is  in  the  set. 

Finally,  if  A  and  B  satisfy  Eq.  (7),  then 

(AB)tK(AB)  >  B^AtKAlB 

-  B^KB  -  K  (9) 


so  that  (AB)  is  in  the  set. 

The  set  of  /f-unitary  matrices  form  a  group,  which  we  shall  call  the 
"Ar-group.  " 

The  If-group  is  infinite  in  the  sense  that  it  contains  an  infinity 
of  elements. 

We  shall  assume  that  the  group  is  parametri zable  at  least  in  the 
neighborhood  of  M  ■  I.  We  can  define  a  set  of  matrices  ,  that  do  not 
depend  either  on  M  or  z,  such  that  we  can  write 

M  =  I  +  Za.Mi  (10) 

where  the  «j’s  are  functions  of  M  which  go  to  zero  continuously  as  M 
goes  to  I.  We  assume  that  this  representation  is  single-valued  for 
sufficiently  small  (X^’s. 

(We  are  skirting  around  some  difficulties,  here.  To  be  proper,  we 
should  show  that  M  is  a  continuous  topological  group,  so  that  we  have 
available  a  proper  definition  of  neighborhood.  In  addition,  the  single¬ 
valuedness  should  be  investigated  further.  However,  we  shall  skip  over 
these  matters . ) 

The  set  are  called  the  "infinitesimal  transformations”  of  the 

group. 
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Consider,  now  a  solution  of  Eq,  (4)  which  reduces  to  the  identity 
at  z  =  z'.  This  will  be  of  the  form  of  Eq,  (10)  with  the  variable 
quantities  being  functions  of  both  z  and  z'— i.e,,  of  position  and  of  the 
position  of  the  boundary  conditions. 


M(z,2')  =  I  +  Sttjlz.z' )Mi 

a.  (2 ,2' )  =  0  at  z  =  z' 


We  will  define  the  quantities 


7i(z') 


B«j(z,z') 

^2 


(11) 


(12) 


The  y.  are  the  coefficients  of  the  first  term  in  the  Taylor  expansion  of 
M(z,z')  with  respect  to  2,  about  the  point  where  M  “  I. 

If  we  substitute  Eq.  (11)  in  Eq..  (4),  we  have 


3M(z,z') 

Bz 


I 


Ba  .  ( 2 ,  z ' ) 

Bz 


M. 


-jRil  +  2a.(z,z')M'} 


We  evaluate  this  at  the  point  z  =  z' ,  when  the  a.  become  zero  and 
find  that 


R(z‘)  =  jly^iz'Mi 

i 


(13) 


Hence R(z)  can  also  be  expressed  in  terms  of  the  set  •  R(z)  does 
not  form  a  group,  but  its  properties  do  depend  on  the  infinitesimal 
transformations  of  the  /(-group  under  discussion. 

We  note  that  these  conclusions  are  equally  valid  whether  or  not  R 
is  constant.  If  the  system  is  nonuniform,  the  coefficients  of  R,  ijy^} , 
are  functions  of  z,  but  otherwise  there  is  no  cliange. 

It  is  a  theorem  of  Lie  groups  that,  given  appiopriate  conditions  of 
single-valuedness,  etc.,  the  infinitesimal  transformations  determine  the 
group  completely. 

From  knowledge  of  the  K-group,  we  will  be  able  to  determine  what 
sets  of  infinitesimal  transformations  we  wish  to  consider.  From  Eq.  (13) 
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we  can  determine  which  of  these  are  actually  involved  in  the  R  describing 
a  particular  system.  If  we  can  so  choose  tlie  complete  set  that  the  subset 
involved  with  the  given  R  are  tlie  infinitesimal  transformations  of  a 
simpler  group,  then  the  system  can  be  studied  in  terms  of  this  simpler 
group . 

C,  THE  ROTATION  GROUP 

We  will  find  that  the  groups  we  are  concerned  with  are  related  to, 
although  not  in  general  identical  with,  the  three-dimensional  rotation 
group.  This  is  perhaps  not  surprising  since  s,  defined  by  Eq,  (.5),  can 
be  regarded  as  a  quantity  analogous  to  the  square  of  the  length  of  the 
vector  X.  That  s  is  invariant  when  X  is  replaced  by  Mx  then  says  that 
the  length  is  constant  under  the  given  operation,  which  may  be  taken  as 
the  definition  of  a  rotation. 

We  call  this  an  analogy  rather  than  a  definition,  because  the  con¬ 
cept  of  length  becomes  obscure.  Since  K  is,  usually,  not  positive 
definite,  s  can  be  zero  or  negative.  Hence  the  “length"  can  be  zero  or 
pure  imaginary,  even  though  the  vector  is  non-null. 

Regardless  of  the  problems  of  geometric  visualization,  it  is  true 
that  these  operators  can  be  regarded  as  rotational. 

We  will,  therefore,  develop  the  properties  of  the  rotation  group 
from  their  abstract  operators. 

The  rotation  group  can  be  described  in  terms  of  the  three  matrices 
which  can  be  taken  as  its  infinitesimal  transformations.  One  possible 
form  for  these  matrices  is 

A, 


Az 


(14) 


(15) 


(16) 
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Note  that  we  are  not  saying  that  Aj ,  Aj ,  and  Aj  are  a  group  —  they 
are  not.  But  they  are  the  infinitesimal  transformations,  i.e.,  the  , 
of  a  continuous  group  of  M  ,  defined  from  them  according  to  Eq.  (10). 

We  will  not  here  be  concerned  with  derivation  of  these  matrices. 

This  can  be  done  by  proper  parametrization  of  the  rotation  group.  We  are 
concerned  with  the  properties  of  the  abstract  group,  not  the  particular 
properties  of  rotation  operators. 

It  can  be  readily  verified  that  these  three  matrices  have  the 
commutation  relations 


[AjAj] 

=  AjAj  “  AjAj  ■  ''2Aj 

(17) 

[A1A3] 

“  “Aj 

(18) 

[A3A3] 

*  Aj 

(19) 

This  is  not  coincidence.  It  is  a  general  theorem  that  the  set  of 
infinitesimal  transformations  of  any  group  have  the  property  that  the 
commutation  of  any  two  of  them  can  be  expressed  in  terms  of  the  set.  By 
forming  linear  combinations  of  the  infinitesimal  transformations-->-which 
is  equivalent  to  a  reparametrization  by  linear  combinations  of  the 
parameters  —  the  commutation  relations  can  be  varied  to  an  extent.  The 
particular  form  of  the  A  in  Eqs .  (14),  (15),  and  (16)  were  carefully 
chosen  to  give  Eqs.  (17),  (18),  and  (19)  which,  as  we  shall  see,  are 
particularly  convenient  relations. 

The  commutation  relations  are  the  fundamental  ones.  We  can,  for 
example,  change  the  basis,  i.e.,  the  coordinate  axes-^of  the  system.  If 
we  do,  the  infinitesimal  transformations  undergo  a  similarity  transforma* 
tion.  They  are  therefore  changed.  But  it  is  easy  to  see  that  the 
commutation  relations  are  invariant  under  such  a  change.  Hence,  the 
commutation  relations  are  independent  of  the  coordinate  system. 

The  commutators  of  Eqs.  (17)  through  (19)  have  remarkable  properties 
in  determining  the  interrelations  of  the  significant  vectors  of  the  system. 

Consider,  first,  an  eigenvector,  of  Aj: 

Vi  =  ^i*i  (20) 
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Then  AjX^,  if  i  *.  is  not  the  null  vector,  is  also  an  eigenvector  of 
Aj  with  eigenvalue  (X.^  +  1).  For,  from  Fq.  (18): 

Aj(A|X^)  =  (AjAj  +  Aj)Xj 

=  (X.  +  l)A,x,  .  (21) 

Hence  we  can  construct  a  sequence  of  eigenvectors  of  Aj  with  eigen¬ 
values  that  differ  by  unity,  that  are  connected  by  Aj . 

Likewise  AjX^ ,  if  it  is  not  the  null  vector,  is  an  eigenvector  of 
Aj ,  this  time  with  eigenvalue  (X.  -  1).  For,  from  Eq.  (19), 

A3(AjXj)  »  (AjAj  -  Ajj)x^ 

=  (X.  -  DAjX.  .  (22) 

Again,  we  can  construct  a  sequence  of  eigenvectors  of  Aj  with  eigenvalues 
that  differ  by  unity,  and  that  are  connected  by  Aj. 

We  still  need  to  show  that  these  two  sequences  are  the  same,  within 
scalar  factors.  Before  doing  this,  however,  consider  the  sequence  formed 
by  Aj.  The  number  of  eigenvectors  of  Aj  is  finite.  Therefore,  the  se¬ 
quence  formed  by  Aj,  starting  with  any  given  eigenvector  of  Aj ,  has  member 


Xjj  with  largest  eigenvalue.  That  is. 

there  must  exist  an  such  that 

A  V  m 

^3*lf 

and 

XjjXjj 

(23) 

^i*if 

0  . 

(24) 

Consider,  now,  the  sequence  formed  from  X^  by  Aj-  Define 

^2^11  "" 

AjXjj.j  “  (25) 
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Then,  using  Eqs.  (25),  (19),  and  (24) 


Likewise, 


*  ^1^2^11 

=  "(AjAj  -  2Aj)Xj^ 


A,x 


11-2 


(AjAj  -  2Aj)*j,_j 

-(2X^jX^  +  2\^_jX^_j) 

~2(\  +  Xjj.j  )Xj,_j 


Let  us  assume  that 


(26) 


(27) 


where 


'  2(A.j,  +  X|,_,  +  ...  \-*  +  i) 


which  we  have  shown  to  be  true  for  1  and  2.  Then 

Aj^Xjj_j_j  = 

=  (AjAj  - 

=  -(AjP^.jXj^,^ j  +  2A.j,_jXj,_j ) 

Hence,  AjX^_j_j  is  proportional  to  ~x^_j  with  scalar  factor 

Ptl-k 


(28) 


(29) 


(30) 


(31) 


Hence,  the  sequences  generated  by  A|  and  Aj  are  identical  within  scalar 
factors. 
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We  should  note  the  possibility  of  more  than  one  sequence.  It  is 
quite  possible  that  there  are  more  than  one.  If  so,  then  the  subspace 
spanned  by  each  set  is  an  invariant  subspace  of  the  whole  space,  and  the 
representation  is  “reducible.”  We  shall  return  to  this  later. 

We  have  used,  so  far,  the  fact  that  there  must  be  a  vector  in  the 
sequence  formed  by  Aj  which  has  the  largest  eigenvalue.  Clearly,  there 
must  also  be  a  smallest  one.  That  is,  there  must  be  a  value  n  such  that 

AjX,  -  (32) 

AjX^  >'  0  .  (33) 

But,  now,  consider 

^1^2*.  =  VjO  =  0  =  -p,.jX.  (34) 

Since  x^  is  not  the  null  vector,  we  must  liave  p^_j  =  0. 

Hence,  the  sequence  of  factors  has  its  two  end  points  fixed. 

We  have,  from  Eq.  (28)  and  the  above,  that 

=  P.-i  =  0  (35) 

with  the  p  connected  by  the  recursion  formula,  Eq.  (31),  where  the  'K  are 
connected  by 

~  ^  '  (36). 

These  conditions  are  sufficient  to  determine  The  solution  to 

Eq.  (31)  which  vanishes  at  k  =  M  is 

+  k{2M  -  2\_j  ) 

+  ’/2  -  MI.2M  -  2^^_i  )  .  (37) 

If  we  let  k  =  m  -  1  and  require  that  P.-j  “  0,  we  can  solve  for  : 

\  =  -T  ■  (38) 

"  2 
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It  is  convenient,  and  conventional,  to  let  the  index  run  over  either 
integral  or  half  integral  values  from  (-l/2)(n  -  1)  to  (+l/2)(n  -  1), 
where  n  is  the  total  number  of  vectors  in  the  sequence.  Then 
M  -  (l/2)(n  -  1)  and  m  =  (“l/2)(n  -  1)  so  the  =  (l/2)(n  “  1)  and 

^  (39) 

+  1)  “  kik  +  1)  (40) 

where 

L  1  1 

“  1)  •••  “1,  0,  +1,  ...  (n  “1)  if  n  is  odd 
^  2 

11  1 

2  ”  1  ~2  ’  ^2  ’  2  ^  even 

Thus,  from  the  commutation  laws  alone,  we  are  able  to  determine  a 
great  deal  of  the  structure  of  the  operators. 

If,  then,  we  can  express  R  in  terms  of  operators  which  obey  these 
commutation  laws,  we  will  from  this  fact,  be  able  to  deduce  many  of  the 
properties  of  the  system.  For  example,  suppose  that  we  can  so  express  R, 
and  suppose  the  boundary  condition,  is  an  eigenvector  of  Aj .  Then  it 
is  evident  that  the  state  vector  must,  for  all  z,  stay  within  the  subspace 
spanned  by  the  set  of  vectors  formed  by  repeatedly  operating  on  with 
Aj  and  Aj  . 

We  must,  therefore,  consider  the  class  of  operators  that  obey  the 
commutation  relations.  These  are  the  '*  representations”  of  the  rotation 
group. 

D.  ANALYSIS  OF  THE  /f-GROUP 

We  have  discussed  in  some  detail  the  three-dimensional  rotation  group 
because  it  is  a  useful  example.  Not  only  does  it  itself  appear,  on  oc¬ 
casion,  but  the  methods  of  analysis  that  are  effective  on  it  can  be  extended 
to  more  general  groups. 

We  do  wish  to  emphasize,  however,  that,  in  the  general  case,  we  cannot 
restrict  ourselves  to  the  three-dimensional  rotation  groups. 
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Consider  an  R  that  is  n  ^  n.  There  are,  then  n*  components  each  of 
which  may  be  complex.  That  R  is  /^-hermitian  restricts  its  form,  but  still 
leaves  n*  degrees  of  freedom.  For  example,  i£  K  =  I  so  that  R  is  hermitian 
we  require  that  the  diagonal  terms  be  real,  and  the  off-diagonal  terms  be 
the  complex  conjugates  of  the  terms  in  the  transposed  position. 

Each  degree  of  freedom  may  vary  with  z  in  any  prescribed  fashion  with¬ 
out  destroying  the  /f-hermitian  property.  Hence  the  most  general  n  ^  n  K- 
hermitian  R(z)  will  involve  n*  arbitrary  functions  of  z. 

To  express  such  an  R  in  terms  of  a  set  of  constant  infinitesimal  trans 
formations  M . ,  as  in  Eq.  (14)  requires  n*  elements  of  Mj  so  that  we  can 
match  the  terms  to  the  prescribed  functions.  The  three-dimensional 

rotation  group,  however,  has  only  three  infinitesimal  transformations.  It 
can,  at  best,  be  used  to  describe  an  R  that  has  only  three  prescribed 
functions  of  z. 

We  must,  therefore,  use  more  complicated  groups  to  describe  systems 
of  greater  complexity. 

There  is  a  useful  qualification  to  this  statement.  Wecanmeasure  each 
mode  of  the  system  against  a  reference  phase  that  varies  in  any  prescribed 
manner  andwith  an  amplitude  scale  factor  that  varies  in  any  prescribed  manner. 
The  transformation  to  such  a  measurement  transforms  R(z),  and  can  be  used 
(at  least  in  principle)  to  remove  2n  of  the  degrees  of  freedom  of  R. 

For  the  general  case,  then,  we  can,  in  principle,  reduce  the  problem 
to  one  involving  (n*  -  2n)  degrees  of  freedom — and  hence,  to  consideration 
of  a  group  that  involves  (n*  -  2rt)  infinitesimal  transformations. 

The  difficulties  of  doing  this  completely  are,  very  often,  formidable. 
For  example,  if  n  =  2,  we  find  we  do  not  need  any  infinitesimal  transfor¬ 
mation.  However,  to  accomplish  the  reduction  we  must  exactly  solve  the 
system.  If  we  could  do  this,  then  the  group  theoretic  analysis  would  be 
trivial  and  unnecessary  anyway. 

It  is,  however,  often  possible  to  accomplish  a  partial  reduction 
without  excessive  difficulty.  For  example,  in  the  2^2  case,  we  again 
have  four  degrees  of  freedom,  initially.  We  can,  then,  represent  the 
2^2  case  in  terms  of  I,  Aj ,  Aj ,  and  Aj ,  where  the  A’s  obey  the  commu¬ 
tation  laws  of  Eq,  (17),  (18),  and  (19).  The  presence  of  the  term  in  I 
prevents  this  from  being  a  representation  of  the  three-dimensional 
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rotation  group.  However,  this  term  is  easily  removed.  In  fact,  we  can 
easily  remove  both  the  I  and  the  Aj  term,  as  we  shall  see  later.  Hence, 
the  problem  can  easily  be  reduced  to  one  involving  the  three-dimensional 
rotation  group. 

Given,  then„ a  representation  of  R  in  terms  of  an  appropriate  set 
of  infinitesimal  transformations  obeying  known  commutation  relations,  it 
is,  then,  possible  to  deduce  much  of  the  structure  of  the  operators  by 
methods  analogous  to  those  of  the  preceding  section.  If,  for  example, 
there  exist  three  operators,  Aj ,  Aj,  and  Aj ,  which  obey  Eqs.  (17),  (18), 

(19),  then  all  the  conclusions  of  the  preceding  section  apply,  even  though 
there  may  be  other  operators  included  in  R  which  do  not  commute  with  A^, 

Aj ,  and  Aj  and  whose  effect,  therefore,  is  interrelated  with  the  effect 
of  the  A  parts  of  R. 

Furthermore,  we  can  separate,  to  an  extent,  the  consequences  of  the 
individual  commutation  laws.  Equation  (18),  for  example,  is  sufficient 
to  couple,  in  chain  fashion,  the  eigenvectors  of  Aj .  That  is,  if  x  is  an  . 
eigenvector  of  Aj ,  with  eigenvalue,  of  X,  then  AjX,  if  it  is  not  the  null 
vector,  is  also  an  eigenvector  of  Aj  with  eigenvalue  (X  +  1).  Aj  is  a 
"raising”  operator,  or,  in  quantum  mechanical  terminology,  is  a  "creation” 
operator. 

Likewise,  if  Eq.  (19)  is  true,  then  Aj  is  a  " lowering”  or  "annihilation 
operator. 

Hence  the  commutation  relations  provide  a  convenient  way  of  obtaining 
the  structure  of  the  system  and  of  describing  the  physical  consequences  of 
the  various  parts  of  R. 

Alternatively,  and  perhaps  more  usefully,  we  can  use  the  desired  struc¬ 
ture  of  the  infinitesimal  transformation  to  find  the  appropriate  representa¬ 
tion  of  R.  If,  for  example,  R  is  obtained  as  the  perturbation  of  a  solved 
Rg ,  then  we  can  use  the  known  eigenvectors  of  Rg  as  the  modes  of  the  system. 
We  can,  then  map  the  couplings  between  these  modes  that  are  introduced  by 
the  perturbation.  We  can,  then,  take  a  chain  of  such  couplings  and  identify 
it  as  Aj ,  and  the  reverse  chain  as  Aj .  There  may  be  several  such  chains, 
each  of  which  is  identified  in  similar  ways.  Each  chain,  then  summarizes 
the  effect  of  one  particular  aspect  of  the  perturbation. 

We  shall,  shortly  discuss  this  process  in  more  detail. 
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Thus,  even  though  the  three-dimensional  rotation  group  is  inadequate 
for  our  purpose,  it  does  serve  as  a  useful  example  and  prototype.  We  shall, 
therefore,  in  the  next  section,  discuss  in  some  detail  its  representation. 

E,  REPRESENTATIONS  OF  THE  ROTATION  GROUP 
1.  General 

The  irreducible  and  non-equivalent  representations  of  the  rotation 
group  are  known  and  may  be  found  in  standard  textbooks.  Our  problem,  here, 
however,  is  to  establish  techniques  so  that  we  may  consider  the  converse 
problem.  That  is,  we  wish  to  write  a  representation  of  the  rotation  group 
in  a  form  that  is  suitable  for  its  identification  with  components  of  a 
given  R. 

The  most  useful  way  of  obtaining  representations  seems  to  be  their 
development  as  dyad  expansions  on  a  set  of  /f-orthogonal ,  maximally  nor¬ 
malized  vectors.  In  the  situation  where  the  given  system  can  be  considered 
as  the  result  of  perturbing  a  system  whose  behavior  is  known-^usually  a 
uniform  system — the  vectors  can  be  the  eigenvectors  of  the  unperturbed 
system.  These  vectors,  then,  are  taken  as  the  "modes”  of  the  system.  The 
perturbation  is  considered  as  introducing  coupling  between  these  modes. 

We  have  introduced  the  concept  of  /f-orthogonality  before,  but  we  will 
review  the  pertinent  conclusions  here. 

We  consider  a  complete  set  of  vectors  (Complete  means  that  they 

span  the  whole  space,  so  that  they  include  a  set  of  linearly  independent 
vectors  in  terms  of  which  any  vector  can  be  expressed.  If  the  space  is 
n-dimensional ,  then  a  complete  set  must  include  exactly  n  vectors,  if  they 
are  all  linearly  independent.  ) 

We  introduce  the  symbology  on  the  indices  (~i )  which  may  be  read  as 
“conjugate  i.”  It  is,  in  other  words,  a  relabelling  of  the  indices.  We 
require  that  this  relabelling  be  “ one - to-one ,  ”  so  that  to  each  vector 
there  corresponds  one  and  only  one  vector  m.^^ ,  which  may  be  the  same  vector, 
and  vice  versa. 

The  set  is  said  to  be  "Af-orthogona  1  ”  if  and  only  if 

-  0  if  i  / 
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(Note:  the  concept  of  conjugate  indices  is  necessary  only  because  K  may 

not  be  positive  definite.  If  K  is  positive  definite,  then  we  must  have 
identical  with  i.  It  is  only  the  indefiniteness  of  K  that  permits 
the  possibility  of  being  different  from  i.) 

The  set  is  said  to  be  "maximally  normalized”  as  well  if 

-  cr. 

where  cr.  =  1  if  i  /  and  is  ±1  if  i  ■  .  In  the  latter  case,  cr^  indi¬ 

cates  the  “parity”  of  the  vector. 

We  have  shown  that,  if  is  f:  coupling,  or  perturbation  parameter, 
such  that  R{<p,z)  is  /(-hermitian  for  all  z,  and  <p,  then  the  eigenvectors 
of  the  uncoupled  system  R(4>,z),  can  be  so  chosen  as  to  be  /(-orthogonal 
and  maximally  normalized. 

Given  such  a  set  of  vectors,  we  can  form  the  dyads  ♦K).  We 

shall  use  these  dyads  to  develop  our  representations. 

Before  proceeding  to  the  more  complex  cases,  let  us  consider,  first, 
a  comparatively  simple  one.  Let 

Aj  =  acr.UjU.^JK  (41) 

Aj  =  Al/oi)cT.u.u^JK  (42) 

where  i  jl  j  so  that  the  vectors  are  not  cross-conjugated  with  each  other, 
and  where  at  is  any  constant,  real  or  complex. 

Then,  from  Eq.  (17),  we  find  that 

^3  =  "  ^;“>“'V>^K)  (43) 

It  is  easy  to  verify,  then,  that  Eqs.  (18)  and  (19)  are  obeyed. 

It  is  of  interest  to  note  that  the  eigenvectors  of  Aj  are  and 
with  eigenvalues  +%  and  respectively.  We  also  note  that 

AjU^  =  0  AjUj  »  -(l/a)er.cr.Uj 

AjU^.  <=  (Xc^.a.u^  XjUj  =  0  (44) 

in  verification  of  our  previous  results. 
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This  representation  is  then  isomorpliic  to  tlie  representation  of  the 
rotation  group  as  2  2  matrices: 


(45) 


Therefore,  if  the  given  R  can  be  expressed  in  terms  of  the  operators 
of  Eqs .  (41),  (42),  and  (43),  then  the  problem  can  be  replaced  by  an  equiv 
alent  one  involving  the  2^2  matrices  of  Rq.  (45). 

As  the  next  simplest  situation,  it  may  be  easily  verified  that  repre¬ 
sentations  of  the  type  of  Eqs.  (41)  to  (43)  can  be  added  if  there  are  no 
cross-coupling  terms.  We  can  take,  for  example, 


Aj  =  a.cr.u.u^/K  + 

Aj  =  -(1/a  .  tR  -  (] /a^  )cr^u^u..^^^K 

We  find  that  these  matrices  satisfy  Eqs.  (17),  (18),  and  (19)  providing 
t,  7,  k ,  and  h  are  all  distinct.  It  is  interesting  to  note,  that  there  is 
no  requirement  that,  for  example  should  not  be  identical  with  or  u..^ 

In  this  case,  there  might  appear  to  be  coupling  between  the  two  sets  of 
vectors,  but  in  fact,  there  is  not. 

We  observe  that  this  representation  is  simply  the  combination  of  two 
representations  of  the  fo'rm  of  Eqs.  (41),  (42),  and  (43). 

The  eigenvectors  of  Aj  are  u.,  ,  n^,  and  ii^.  If  the  R  of  a  given 

system  can  be  expressed  in  terms  of  Aj,  A2,  and  A3,  then  the  subspace 
spanned  by  u^.  and  u^.  ,  and  tliat  spanned  liy  and  ,  are  invariant  sub¬ 
spaces  of  the  system.  Eacli  subspace  can,  then,  be  replaced  by  an  anal¬ 
ogous  2  2  problem  involving  the  A’s  of  Eq.  (45). 
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Clearly,  we  can  add  as  many  representations  together  as  we  wish, 
providing  they  are  wholly  decoupled. 

Let  us  consider,  now  a  representation  of  higher  order,  We  wish  Aj 
to  be  an  operator  which  establishes  a  sequence  of  vectors.  Consider  for 
example : 

Aj  -  +  0i2cr^.U^.U.^j^K 

where  i  /  j  /  k.  Then 

AiUj  »  0 

AjU.  - 

A,u^  =  a^cr.ai^u. 

so  that  Aj  establishes  the  sequence  k  j  i . 

We  wish  A2  to  establish  the  same  sequence  in  reverse  order.  We  there¬ 
fore  let 

so  that 

AjU.  -  yio-^cr.u. 

AjU^.  -  72‘^yCr*«i 

AjUj  «=  0  .  (48) 

We  still  need  to  determine  the  relationship  between  the  c '.istants.  From 
Eq.  (17)  we  find 

-  aj7jO-.u.u.^.  ♦K  -  a2r2‘^y“y'''v;  .  (49) 

We  can  easily  find  that 

[AjA,]  -  ailaiTi  -Y®272)^.“<“~/*K  +  «2(®272  -  . 

(50) 
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For 

this 

to  e<iual 

“Aj  according  to  Eq 

.  (18),  we  must  have 

-1 

= 

-1  • 

(51) 

This  requires 

that 

-2 

Hence  we 

take 

(52) 

=  -(2/aj  )cr^.u^.u,^^  ♦K 

-  {2/a, 

(53) 

Aa 

=  a.u.u,,.tK  -  cT^u^ 

'K  . 

(54) 

It  can  be  easily  confirmed  that  Eq.  (19)  is  satisfied. 

The  eigenvectors  of  Aj  are  U.,  U^,  and  with  eigenvalues  +1,  0, 

"1,  as  expected. 

This  representation  is  isomorphic  to  the  3^3  representation  of  the 
rotation  group,  for  which  we  may  use  the  matrices  given  in  Eqs.  (14), 

(15),  and  (16). 

Clearly,  again,  we  can  add  representations  together  providing  they 
do  not  have  vectors  in  common. 

The  same  technique  can  be  used  to  develop  representation  coupling 
any  number  of  vectors  together.  We  take  Aj  as  the  sum  of  scalars  times 
the  dyads  that  move  us  along  the  chosen  sequence.  We  take  Aj  as  the  sum 
of  different  scalars  times  the  dyads  that  move  us  along  tlie  same  sequence 
in  the  reverse  order.  We  can,  then,  compute  Aj  in  general  terms  from 
Kq.  (17).  From  either  Eq.  (17),  or  (18)  we  obtain  the  necessary  relations 
between  tlie  scalars  of  Aj  and  those  of  Aj .  This,  then,  fixes  Aj  ,  Aj  ,  and 
Aj  ,  within  the  allowable  variations,  so  that  Eqs.  (17),  (18),  and  (19) 
are  satisfied.  We  have,  then,  a  suitable  representation  of  the  rotation 
group . 
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In  this  way,  we  are  able  to  obtain  a  variety  of  representations  of 
the  rotation  group,  and  to  discuss  them  in  terms  of  the  coupling  of 
vectors . 

2.  Example 

To  illustrate  our  discussion  of  the  meaning  of  the  different  repre¬ 
sentations,  we  shall  consider  a  simple  case. 

Let  us  consider  the  system  matrix  of  a  transmission  line: 


/  0  jijZ 

\^IZ  0 


for  which  a  suitable  K  is 


K 


(55) 


(56) 


We  are  not  here  assuming  that  Z  is  independent  of  z  so  that  this  may, 
itself,  be  a  non-trivial  problem.  We  shall,  however,  show  the  relation 
between  this  problem  and  one  of  higher  dimensionality,  and  therefore 
greater  apparent  complexity. 


If  we  take 


(57) 


(58) 
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Hence 


R 


(59) 


What  we  want  to  do  now  is  to  consider  whether  or  not  there  are 
situations  of  greater  apparent  complexity  which  are,  nevertheless,  es¬ 
sentially  the  same  as  Eq.  (55),  and  for  which  a  solution  to  Eq.  (55)  can 
be  used.  This  will  be  true  if  Eq.  (59)  is  true,  but  with  an  Aj  and  Aj 
that  is  more  complex  than  Eq.  (58). 

To  consider  this  problem,  let  us  proceed  in  the  opposite  direction 
by  considering  a  more  complex  vector  constructed  from  the  vector  of 
Eq.  (55).  Let  us  consider: 

y 


bXjXj 

cxl 


(60) 


where  a,  b,  and  c  are  constants  to  be  determined.  This  form  is  a  general¬ 
ization  of  the  “direct  product"  or  Kronecker  product  of  y  with  itself. 

We  see,  then,  that 

/  2aXj (d*j/dz)  \ 

=  (  6*,(dXj/dz)  +  6xj(dx,/dz) j 

\  2cXj(dXj/dz)  / 


/  ?.ax^{-jfiZx^) 

f  bx^(-j(3x^/Z)  +  bx^{-jfiZx^) 
\+  2cx^(-jfix^/Z) 


(0  2a^Z/b 

bp/aZ  0 

0  2c/S/cZ 


The  conservation  law  states  that 

S  »  **K*  =  ^ 


(61) 


(r.2) 


SO 


is  invariant.  We  require  the  existence  of  o  conservation  law  for  the  y 
vector.  We  cajt  obtain  such  a  law  by  observing  that 


5 


2 


+  2x, 


(63) 


must  also  be  invariant. 

For  this  to  be  a  conservation  law  for  the  y  vector  we  require, 
within  a  constant  factor,  that 


a  ■=  ei'f'f/A* 

b  =  A 

c  =  . 


(64) 


Then,  if  we  set 


and 


R' 


K' 


0  0/1 
0  1  0 

4*  0  0 


0 

/24  '^/8/7  0 

0  /iti't’pfz/k  0 


we  can  easily  see  that  R*  is  /f* -hermitian. 
We  can,  now, set 


a; 


^1 


'0  ot,  0 

0  0  a, 

^0  0  0 


(65) 


(66) 


A' 


a; 


1  0 
0  0 


0 

■1/a, 

0 

0 
0 


0 

0 

"1/a, 


0 

0 

0 


(67) 


0  0  -j 


and  find  that  these  matrices  obey  the  conservation  rules. 
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If  we  set 


a 


1 


a 


(68) 


then  we  find  that 


R' 


K 


(69) 


which  is  the  same  form  as  Eq.  (59). 

The  point  of  this  development  is  not  the  determination  of  R'  from  R, 
as  we  have  done,  but  rather  the  converse. 

Given  a  three-dimensional  system  with  the  matrix  R'  of  the  form  of 
Eq.  (66),  where  Z  and  /S  may  be  functions  of  z.  We  do  not  need  to  solve 
it  directly.  Instead,  it  is  sufficient  to  solve  the  equivalent  two- 
dimensional  system  with  the  R  of  Eq.  (55).  The  solution  of  the  three- 
dimensional  system  is,  then,  the  vector  y  which  is  determined  through 
Eq.  (60),  from  X,  the  solution  of  the  two-dimensional  problem. 

This  technique,  titen  offers  hope  that  a  multidimensional  problem  can  be 
reduced  to  one  or  a  sequence  of  problems  of  lower  dimensionality  and  the 
answer  obtained  by  the  appropriate  direct  products  of  the  solutions  of 
the  simpler  prol)lems. 

F.  CONSTANT  R 

Since  a  uniform  system  (i.e.,  one  such  that  R  is  of  simple  structure 
and  constant)  is,  in  principle,  always  soluble  by  the  eigenvectors  (and, 
if  necessary,  the  generalized  eigenvectors)  of  R,  the  group  theoretic 
properties  of  such  systems  are  principally  of  academic  interest. 
Consideration  of  such  systems,  however,  may  clarify  the  significance  of 
some  of  the  statements  that  we  have  made. 

There  are  several  ways  of  stating  the  property  that  makes  uniform 
systems  directly  soluble.  We  may  observe,  first,  that  the  eigenvectors 
themselves  each  form  a  one-di mens iona I  representation  of  the  rotation 
group.  This  is  the  trivial  representation  of  the  rotation  group  that 
assigns  to  each  operator  the  scalar  ]. 
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In  more  conventional  terms,  we  know  that  we  can  use  the  eigenvectors 
as  a  basis  so  that  R  becomes  diagonal.  The  appropriate  similarity  trans¬ 
formation  is  a  constant  one,  and  the  system  has,  then,  been  resolved  into 
/1-orthogonal  and  maximally  normalized  modes.  It  is  these  modes  that 
permit  direct  solution. 

A  third  way  of  putting  it  is  to  say  that  the  eigenvectors  each  gen¬ 
erates  a  one-dimensional  subspace  that  is  invariant  under  R  and  therefore 
under  any  M  that  is  a  solution  of  the  system. 

When  we  look  for  two-dimensional  representations,  we  are  in  fact 
looking  for  two-dimensional  subspaces  that  are  invariant  under  R.  Then, 
if  X  initially  lies  in  such  a  subspace,  it  continuously  transforms  into 
a  vector  in  the  subspace,  and  hence  is  retained  permanently  within  it. 

The  subspace,  again,  is  invariant  under  the  system  matrix,  M. 

A  change  to  such  a  basis  transforms  R  into  a  quasi -diagonal  form  with 
a  2  2  matrix  on  the  diagonal  in  the  position  corresponding  to  the  given 

subspace.  If  the  basis  is  such  that  the  entire  space  is  divided  into  dis¬ 
joint  invariant  subspaces,  the  transformed  R  is  quas i -di agonal  with  2^2 
matrices  along  the  diagonal. 

Any  subspace  that  is  spanned  by  any  two  eigenvectors  of  R  is  evi¬ 
dently,  an  invariant  subspace.  On  the  other  hand,  if  the  eigenvectors  are 
not  degenerate,  such  a  subspace  can  be  identified  by  two  vectors  that  are 
not  eigenvectors  of  R. 

It  is  not  evident  that  any  two-dimensional  invariant  subspace  is 
necessarily  a  representation  of  the  rotation  group.  One  finds,  generally, 
an  appropriate  determination  of  Aj  and  Aj.  However,  there  is  ho  guarantee 
that  the  resultant  Aj,  determined  from  Eq.  (17),  will  "fit”  the  remaining 
terms.  If  not,  however,  we  can  correct  the  discrepancy. 

It  is  shown  in  the  appendix  that  we  can  make  a  transformation  of  x, 
and  hence  of  R,  that  eliminates  all  " self-coupling”  terms  —  i.e,,  all  terms 
involving  Hence  we  can,  after  such  a  transformation,  always 

write  any  two-dimensional  subspace  in  terms  of  Aj  and  Aj  alone. 

The  transformations  that  are  involved  amount  to  a  z-dependent  renor¬ 
malization  of  the  modes  of  the  system.  .Specifically,  we  consider  the 
system  that  would  occur  if  there  were  no  c ross -coupl i ng  terms  —  i.e.,  if 
all  terms  involved  cr  ,u  ^K(  i  /  j)  were  eliminated.  The  solution  to 
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this  systorn  is  iiscmJ  as  the  reforonre  for  tho  full  -system,  the  phase  angles 
and  amplitudes  of  the  modes  being  measured  by  reference  to  the  solution 
of  the  non -c ross -coup  led  system. 

This  transformation  is  valid  not  only  for  the  uniform  systems  con¬ 
sidered  in  this  section,  but  also  for  non-uniform  systems.  It  is  always 
possible  to  eliminate  completely  all  se 1 f -coupl i ng  terms,  regardless  of 
their  z -dependence . 

G.  NON-UNirOHM  SYSTKMS 

If  R  is  a  function  of  z,  then  we  have  a  different  situation. 

If  the  z-dependence  is  only  in  the  eigenvalues  of  R,  then  the  eigen¬ 
vectors  can  be  chosen  to  be  constant,  and  the  previous  considerations 
apply  without  exception. 

If,  however,  it  is  not  possible  to  choose  eigenvectors  that  are 
constant,  then  we  can  no  longer  use  the  one-dimensional  representations. 

In  general,  in  this  case,  we  must  operate  with  invariant  subspaces 
of  R  that  are  independent  of  z.  If  we  attempt  to  use  subspaces  that  are 
z-dependent,  we  can  still  find  operators  with  the  proper  commutation 
rules,  in  terms  of  which  R  can  be  decomposed.  However,  these  operators 
will  not  be  independent  of  z,  and  it  does  not  generally  follow,  then, 
that  M  can  also  be  expressed  in  these  terms.  It  will  follow  only  if  the 
derivatives  of  ,  and  Aj  are  also  decomposable  in  terms  of  Aj ,  Aj , 

and  Aj.  This  appears  to  be  a  rather  stringent  condition,  although  it 
may  lead  to  some  special  situations  of  interest. 

That  constant,  invariant  subspaces  exist  is  evident.  At  worst  we 
can  take  the  whole  space,  which  is  constant. 

Of  probably  greater  i uteres t  a  re  those  cases  in  which  there  exist 
proper  invariant  subspace.s.  Assuming  it  i  .s  possible  to  express  R,  pos¬ 
sibly  after  adjiistment  of  the  reference  phases,  in  terms  of  constant  A’ s 
within  the  subspae.e,  with  the  z-dependenre  rontnined  wholly  within  the 
y  factors,  then  R  will  bo  i.somorpliic  two-d  i  mens  i  onn  I  representations  of 
the  rotation  group.  The  problem,  then,  can  be  replaced  by  a  simpler  one. 
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H.  DECOMPOSITION  OF  dn/dz  INTO  A  COMMUTATOR 


As  an  immediate  application  of  some  of  these  concepts,  we  will  con* 
aider  the  case  when  R,  which  is  not  necessarily  constant  but  has  the  property 
that  there  exists  a  constant  matrix  C  such  that 


dz 


C«  -  BC 


[C.  R] 


(70) 


It  was  shown  in  the  Second  Interim  Report  that  this  class  of  non- 
uniform  systems  has  the  explicit  solution 


M(z)  «  expI-jR  “  Oz  exp  Cz 


(71) 


A  problem  that  was  unanswered  at  that  time  was  the  determination  of 
the  systems  that  belong  to  this  class,  what  their  general  properties  are, 
and  what  tests  can  be  made  to  determine  if  a  given  system  belongs  to  the 
class.  Also,  no  systematic  procedure  was  available  for  the  determination 
of  the  constant  C  matrix. 

We  shall  find  that  the  representation  of  R  in  terms  of  the  infinites¬ 
imal  transformations  Aj ,  Aj ,  and  Aj  of  the  rotation  group,  when  this  is 
possible,  allows  us  to  find  an  answer  to  these  questions  in  certain 
circumstances. 

We  shall  here  confine  our  attention  to  the  case  when  the  system  can 
be  expressed  as  a  trajectory  in  a  representation  of  the  three-dimensional 
rotation  group.  That  is,  we  will  assume  that  there  exists  constant  oper¬ 
ators  Aj ,  Aj,  and  Aj  such  that  Eqs.  (17)  through  (19)  are  obeyed  and  such 
that  R  can  be  written  as 


R 

so  that 

dz 

It  appears  likely  that  other,  more  complicated  situations  could  be 
handled  in  similar  fashion. 


7jAj  +  72^2  ^  T'jAj 


dK ,  dk  dk . 


(72) 
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It  is  reasonable,  tlien,  to  look  for  a  constant  matrix  C  which  can 
be  represented  in  the  same  way; 


C  =  01, Aj  +  oijAj  +  oijAj  .  (73) 

To  C  we  can  add  any  other  matrix  that  commutes  with  R,  or  that  commute 
separately  with  A,,  Aj ,  and  Aj . 

The  commutator  then  is 

[C,  R]  »  Xotj-yj  -  0(^7, )  [A, ,  Aj] 

+  j(o(,73  -  0(37, )[Ai,  A3] 

+  j  (0(273  -  ^3^ 

=  -2;(o(,72  -  "273)^3 
-j(o(i73  - 

+J(a273  “  “372)^2  .  (74) 

Comparing  Eqs.  (72)  and  (74),  we  see  that  the  necessary  and  suffi¬ 
cient  conditions,  under  the  given  assumptions,  for  the  existence  of  a 
decomposition  into  a  commutator  are  that 


d7, 

dz 


«i73  +  «3ri 


dy^ 

—  -  a273  -  «372 


dy^ 


(75) 


(76) 
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If  tlie  y  are  solutions  of  Kq.  (76),  atul  then  Kq.  (70)  holds.  This 
does  not  def)end  on  the  constancy  of  the  (X,  and  therefore  is  the  general 
condition  for  Kq.  (70). 

If  the  a  are  required  to  be  constant,  then  Eq.  (76)  has  the  solution 


r(z)  =  exp(Pz)r(0) 


(77) 


where 

'■  ■(?) 

and  P  is  the  matrix  of  Kq.  (76). 

We  can  find,  then,  by  the  usual  methods  that 

2Q!j0t2  +  otjX  sin  Kz  \  2o(J(l  -  cos  xz) 

+(2o(ja2  -  K^)  cos  Kz / 

-  cos  Kz)  /  2nr,jn'.2  -  a^K  sin  kz 

\  +(2a,a2  -  K*)  cos  KZ 


eP‘  -  -2 


=  K‘ 


3rt2K  sin  Kz 


"2(XjK  sin  KZ 


K*  - 


+2«jaj(l  -  cos  kz)/  \  +2aj«j(l  -  cos  kz)/  \  +4rtjaj  cos  kz 


(78) 


(79) 


where 


^«l«2  ■  «3 


and  may  be  real,  imaginary,  or  complex. 

This  is  valid  providing  k  /  0.  If  k  =  0,  the  solution  is 


/  (1  +  ajz/2)* 


i.2 


atz 


-a,z(l  +  a^z/2)\ 

(1  -  ajz/2)^  «2z(l  -  cc^z/2.) 

\  2ajZ(l  +  ajz/2)  -2a,z(i  -  ajz/2)  1  -  o»zV2  / 


alz^ 


(80) 


(81) 
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These  solutions  are  easily  obtained  by  considering  the  power  series 
expansion  of  the  exponential  and  using  the  Cay  ley -Hami 1  ton  Theorem  which, 
in  this  case,  tells  us  that 


p3  „  _^2p 

Therefore,  a  necessary  and  sufficient  condition  under  the  given  assumption 
is  that  the  y  vary  according  to  Eq .  (77)  where  exp{Pz)  has  the  form  either 
of  Eq .  (79)  or  (81).  Such  a  set  of  y  generates  a  class  of  systems  that 
can  be  solved  by  the  method  used  here.  Furthermore,  the  matching  to 
Eq.  (79)  or  (81)  will  determine  the  a,  and  therefore  the  A. 

As  a  test  for  whether  a  given  system  belongs  to  this  class,  the 
formulas  given  leave  something  to  be  desired.  We  can  develop  some  tests 
which  are  necessary,  although  not  sufficient,  and  which  are  far  easier 
to  apply. 

We  can  consider  Eq.  (75)  as  determining  the  (X  in  terms  of  the  y. 

We  observe  that  the  determinant  of  the  right-hand  side  vanishes.  There¬ 
fore  these  equations  must  be  redundant.  We  can  find  that  the  condition 
for  this  is  that 


271  • 


<^72 

dt 


<^7, 


^3- 


di 


+  272- 


d7, 

dt 


(82) 


or  that 


47,72  “  7l  +  c 


(83) 


where  c  is  a  constant. 

Equation  (83)  is  a  necessary  and  sufficient  condition,  if  R  and  C 
can  be  represented  in  terms  of  the  rotation  group,  for  dR/dz  to  be  ex¬ 
pressible  as  a  commutator  at  all,  although  it  is  not  sufficient  to  as¬ 
sure  the  constancy  of  C.  We  presume  that  Eq.  (83)  is,  in  some  way, 
equivalent  to  the  condition  that  the  eigenvalues  of  R  be  constant. 

We  can  also  eliminate  from  Eqs.  (75)  any  two  of  the  y.  We  find 
then,  assuming  that  the  (X  are  constant, 

{d^/dz^  -  a^d/dz  -  ia^a^a^)y^  =  0  .  (84) 
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It  is  ner  pssary,  a  I  Lliougli  not  sufficient,  that  each  of  the  y  shall 
obey  an  e(|uation  of  the  form  of  K(|.  (84)  with  the  same  coefficients  for 
a  constant  C  to  exist. 
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II  MATHEMATICAL  TECIINIOUES  FOR  NON-UNIFORM 
COUPLED- MODE  SYSTEMS 


A.  INTRODUCTION 

In  Part  II-A  of  Interim  Report  2,  we  considered  systems  described 
by  the  vector  differential  equation 


dx(z) 

dz 


-jR(z)xiz) 


(85) 


where  R  may  be  a  function  of  z.  We  considered,  in  particular,  those 
systems  in  which  R  is  related  to  a  matrix  A  by  the  equation 


dz 


-jtA,  R] 


->(AR  -  RA) 


(86) 


s<{uare  brackets  being  used  for  the  commutator. 

Interest  in  such  systems  was  stimulated  by  the  fact  that,  if  A  is 
constant — i.e.,  independent  of  z  —  then  tlie  system  is  soluble.  Thus,  this 
is  a  class  of  non-uniform  systems  for  which  it  is  possible  to  obtain  an 
exact  and  explicit  solution. 

Our  purpose  here  is  to  study  this  situation  in  greater  depth  so  as 
to  obtain  better  understanding  of  what  is  going  on,  and  to  determine,  if 
possible,  how  to  generalize  the  class  of  soluble  systems  in  some  useful 
manner. 

We  may  note,  at  the  start,  an  interesting  and  suggestive  symmetry. 
.Suppose  there  exists  a  matrix,  K,  which  we  call  a  "metric”  and  which  is 
constant,  non-singular,  and  hermitian  (K  =  K^,  the  dagger  meaning  the 
complex  conjugate  transpose).  And  suppose  that  R  is,  at  all  values  of  z, 
^-hermitian  so  that 

KR  =  R<K  •  (87) 

[The  general  significance  of  this  relation  we  have  discussed  elsewhere. 
The  pertinent  aspects,  will  be  reviewed  later.  For  the  moment,  let  us 
simply  take  Kq.  (87)  as  given.] 
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Let  us  define  the  matrix  X  as 


X  =  m*K 


(88) 


Then,  since  K  is  constant 


dz 


— 

dz 


+ 


=  -ylUx^K  +  jx*^R*K 


=  -jRx**K  +  ;xx»im 


«  -jRX  +  ;XR 

“  -jtR,  X]  (89) 


Equation  (89)  is  of  the  same  form  as  Eq.  (86),  but  with  X  in  place  of  R 
and  R  in  place  of  A.  We  may  suspect  then,  that  Eq.  (86)  is  a  useful  al¬ 
ternative  representation  of  the  system  represented  by  Eq.  (89).  If  A  is 
constant,  it  is  certainly  simpler. 

Even  if  A  cannot  be  chosen  as  constant,  it  may  be  possible  to  so 
choose  A  that  Eq.  (86)  is  still  a  considerable  simplification. 

This  line  of  thought  suggests  the  possibility  of  continuing  the 
process.  If  A  cannot  be  chosen  as  constant,  is  it  possible  to  choose  it 
so  that 


dz 


-y[B.  A] 


where  R  is  constant?  Or  if  not,  can  we  choose  R  so  that 


(90) 


dB 

dz 


->tC,  R] 


(91) 


where  C  is  constant?  And  so  on. 

If,  at  any  point  in  this  process,  we  do  obtain  a  constant  matrix, 
then  we  can  start  from  that  point  and  try  to  work  our  way  back  up  the  line. 
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If,  for  example,  C  is  constant  in  Kq.  (91),  then  the  solution  to  Eq.  (91) 
i  s 

B  =  .  ('’2) 

This  expression  for  B  can  then  be  substituted  in  Eq.  (90)  and  a  solution 
for  A  sought.  And  so  on. 

This  suggests  that  we  have,  here,  a  powerful  tool  for  the  analysis 
of  non-uniform  systems.  Furthermore,  it  provides  a  classification  of  such 
systems  according  to  the  number  of  steps  involved  before  a  constant  matrix 
is  reached. 

The  results  that  will  be  reported  here  do  not  completely  establish 
this  method  of  analysis.  There  remain  several  important  questions  for 
which  we  have  been  unable,  so  far,  to  find  rigorous  answers.  In  partic- ■ 
ular,  we  may  cite  the  following: 

(1)  Given  a  matrix  R  which  can  be  expressed  as  Eq.  (86).  I® 
always  possible  to  choose  A  so  as  to  obtain  Eq.  (90) — i.e., 
which  has  constant  eigenvalues?  It  appears  that  we  always 
can,  at  least  for  R  K-hermitian.  But  we  have  not  proven  it. 

(2)  Assuming  that  the  answer  to  the  first  question  is  yes,  does  the 
sequence  of  Eqs.  (86),  (90),  (91)  always  lead  eventually  to  a 
constant  matrix,  or  can  it  be  made  to?  We  suspect  that  it  does 
not  —  that  it  is  possible  to  have  an  R  such  that  the  sequence  of 
equations  cannot  be  made  to  terminate.  However,  this  is  pure 
speculation. 

(3)  Assuming  that  the  answer  to  the  second  question  is  no,  then  in 
what  sense  is  Eq.  (90)  simpler  than  Eq.  (86),  Eq.  (91)  simpler 
than  Eq.  (90),  etc.  Is  it  a  valid  approximation  scheme,  for 
example,  to  take  R  >  R^  as  the  first  approximation,  A  =  Ag  as 
the  second,  etc?  We  suspect  that  it  is,  but  we  have  not  proven 
convergence,  or  studied  the  rate  of  convergence  to  determine 
its  usefulness. 

There  are,  then,  serious  questions  remaining.  Nevertheless,  the  results 
we  have  obtained  are  sufficient  to  indicate  the  value  of  the  approach. 

B.  A-CONSERVATIVE  SYSTEMS 

Before  considering  the  problem  itself,  we  shall  review  some  of  the 
pertinent  implications  of  the  existence  of  a  K  metric. 
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We  restrict  our  attention  to  systems  that  exhibit  a  quadratic  con¬ 
servation  law.  That  is,  we  assume  that  there  exists  a  constant  non¬ 
singular  matrix  K,  which  can  be  taken  as  hermitian,  such  that  the  scalar 
quantity,  s, 

s  =  x*Kx  (93) 

is  preserved  (i.e.,  is  independent  of  z)  if  the  vector  is  any  solution  of 
the  system. 

The  conservation  of  energy,  the  linearized  Manley-Rowe  relations,  and 
Chu’s  kinetic  power  theorem  are  examples  of  such  a  conservation  law  and 
lead  to  such  a  K.  It  appears  that  any  system  t}>at  does  not  contain  loss 
mechanisms  has  such  a  law,  although  this  has  not  been  generally  proven. 

At  least  it  can  be  said  that  the  class  of  systems  with  such  a  law  does 
include  a  great  many  systems  of  considerable  interest. 

If  we  differentiate  Eq.  (93)  and  substitute  Eq.  (8.“)),  we  see  that  a 
sufficient  condition  is  that  R  be  /(-hermitian,  or  that  Eq.  (87)  hold. 

The  necessity  of  this  requires  a  little  more  subtle  analysis,  but  can  be 
shown,  also. 

We  have  also  shown  that  Eq.  (87)  implies  that  the  eigenvectors  and, 
if  necessary,  the  generalized  eigenvectors  of  R  can  be  chosen  so  that  they 
are  "/(-orthogonal  and  maximally  normalized.”  That  is,  if  the  set  of 

vectors  {u j }  are  the  eigenvectors  and  generalized  eigenvectors  and,  if 
any  degeneracies  among  them  are  suitably  resolved,  then  they  can  be  so 
normalized  that 


=  cr.8..  (94) 

where  cr  ^  =  ±1  and  S.j  is  the  Kronecker  delta  and  i  s  1  i  f  i  •=  j,  or  0 
otherwise.  Tlie  .symbol  (~j )  which  can  be  read  "conjugate  j”  indicates  a 
possible  relabelling  of  the  indices.  For  a  given  u^,  tliere  is  only  a 
single  vector  of  the  set  such  that  the  scalar,  Eq.  (94),  does  not  vanisli. 

We  call  this  vector  .  It  may  be  the  same  as  u^,  or  it  may  be  different. 

The  complication  of  the  conjugate  indices,  and  the  parity  terms,  cr^  , 
arise  because  K  may  not  be,  and  in  general  is  not,  positive  definite. 

In  the  analysis  that  follows,  we  will  assume  that  R  is  K-hermitian 
and  that  we  know,  or  can  find,  a  complete  set  of  vectors  {Uj(z)}  that  are. 
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at  all  2,  eigenvectors  and  generalized  eigenvectors  of  R,  and  that  are 
everywhere  /f-orthogona  1  and  maximally  normalized,  and  such  that  the  crj 
are  everywhere  constant. 

We  shall  also  assume  that  these  vectors  can  be  differentiated  as 
many  times  as  necessary.  This,  in  effect,  assumes  not  only  that  R  can 
be  differentiated  but  that  its  “structure”  remains  constant.  If,  for 
example,  R  includes  a  description  of  the  transverse  modes  of  an  electron 
beam,  we  must  require  that  the  longitudinal  field  shall  not  reverse  sign, 
since  this  would  reverse  the  parities  of  the  transverse  modes. 

C.  GENERAL  THEOREMS 

We  shall,  now,  prove  some  theorems  that  are  of  value  in  establishing, 
partially,  this  method  of  analysis,  and  that  help  to  give  insight  into 
its  significance. 

Theorem  1  —  A  necessary  and  sufficient  condition  for  the  existence  of 
a  matrix  V,  constant  or  not,  that  satisfies 

dU 

— -=  ->tV,  U]  (95) 

dz 

is  that  the  eigenvalues  of  U  be  constant. 

This  theorem  was  proven  in  Part  I1*A  of  Interim  Report  2,  and  we 
shall  not  repeat  its  proof  here. 

If  we  consider  U  •=  R,  this  appears  to  be  a  restriction  on  the  system 
being  considered.  It  may  not  be,  however.  If,  in  Eq.  (85),  we  let 

X  «  Sy 

where  S  is  a  function  of  z,  then 

f/y  ,  ,  (/S 

— ^  =  -j  S-’RS  -  y  (96) 

(iz  dz 

Hiul  tlie  rjiicstion  now  is  whether  we  c^n  so  choose  S  that 

IQ 

R'  =  S-’RS  -  >8"'-^-  (97) 
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has  constant  eigenvalues.  We  suspect  that  this  is  always  possible,  at 
least  in  principle.  If  so,  then  the  entire  class  of  /^'Conservative  non- 
uniform  systems  can  be  studied.  This  is,  however,  only  a  conjecture  at 
the  present  time. 

We  now  seek  a  relation  between  the  V  of  Eq.  (95)  and  the  detailed 
structure  of  U.  It  is  provided  by  the  following  theorem: 

Theorem  2  —  If  we  have  a  complete  set  of  {Uj}  of  eigenvectors  and 
generalized  eigenvectors  of  U  such  that 

dUj 

— -  -  ->Vu,  (98) 

dz 

then  Eq.  (95)  is  satisfied  with  this  matrix  V.  Conversely,  given  Eq.  (95), 
we  can  choose  a  complete  set  of  eigenvectors  and  generalized  eigenvectors 
of  U  that  satisfy  Eq.  (98). 

To  prove  the  first  part,  first  consider,  an  eigenvector  of  1): 

Ihij  -  XjUj  (99) 

By  Theorem  1,  is  constant.  Hence,  differentiating,  we  obtain 


or 


+  V 


dz 


K 


dz 


mu,  -  jYh.m, 


-  jUVuj  -  yvuoj 

=  -j[V,  U]u.  .  (100) 

Consider,  now,  a  generalized  eigenvector  such  that 

Un.  «  \jn.  +  n._j  (101) 
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Again  differentiating 


dU 

dz 


du .  du 

-u—  +  A; 


dz 


dz 


dz 


=  >UVu.  -  jV(A,u,  +  u..,) 

=  jUVu.  -  jVUu. 

•=  -jtV,  U]u^  .  (102) 

Hence,  Eq.  (100)  holds  for  any  u^,  eigenvector  or  generalized  eigenvector. 
Since  it  holds  for  a  complete  set  of  eigenvectors,  the  operators  must  be 
equal  and  Eq .  (95)  must  follow. 

To  prove  the  second  part,  we  again  consider  first  an  eigenvector  and 
differentiate  Eq.  (99),  substituting  now  from  Eq.  (95).  We  find 


(U  -  ^,1) 


+ 


0 


(103) 


Alternatively,  if  we  consider  a  generalized  eigenvector  and  dif¬ 
ferentiate  Eq.  (101),  we  find 


(U  -  ^J) 


(/ii 

t-i 

dz 


+ 


jVu;., 


(104) 


Equation  (103)  is  an  eigenvector  relation  and  Eq.  (104)  is  a  gen¬ 
eralized  eigenvector  relation.  That  is,  for  example,  Eq.  (103)  states 
that  the  vector  (du^/dz  +  JYu .)  must  be  in  the  subspace  spanned  by  the 
eigenvectors  of  U  with  eigenvalues  k.. 

If  there  is  degeneracy,  we  can  so  resolve  the  degeneracy  that 
{du./.dz  +  jVu  ^ )  is  a  scalar  function  times 


dz 


whether  is  an  eigenvector  or  a  generalized  eigenvector. 


(105) 
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Let  lis  now  define  a  new  set  of  vectors  hy  the  relation 


u.  ■=  e  V.  . 


If  we  let 


=  -j 


(106) 


(107) 


then  the  right  si<Je  of  Lq.  (10!))  vsinishes  and  we  are  left  with 


tlz 


-;Vv. 


(108) 


Hence,  Kq.  (OS)  does  imply  that  it  is  possible  to  choose  the  set  of 
eigenvectors  and  generalized  eigenvectors  {u^}  so  that  Kq.  (Ofl)  is  true. 

We  note  that  we  have  a  considerable  variation  allowed  to  us,  both 
in  the  set  {u^}  and  in  the  corresponding  matrix  V.  As  regards  the  set 
of  vectors,  any  degeneracies  may  be  resolved  arbitrarily,  and  the  manner 
of  resolution  may  be  itself  a  function  of  z.  Each  eigenvector  may  be 
multiplied  by  an  arbitrary  non>vani shi ng  scalar  function  of  z  without  af¬ 
fecting  its  being  an  eigenvector.  But  to  each  such  set  of  vectors  there 
is  a  particular  V. 

Conversely,  we  can,  in  Eq .  (9.5)  add  to  V  any  matrix  which  commutes 
with  U  without  essentially  changing  Eq.  (95).  In  particular,  we  can  add 

any  suitable  function  of  U  [a  function  of  f{x,  z)  is  suitable  if  for  all 
z,  it  is  analytic  with  respect  to  *  in  a  region  that  includes  all  the 
eigenvalues  of  U.  We  may  then  add  to  V  the  matrix  /(U,  z).]  The  corres¬ 
ponding  set  of  vectors  is,  then,  multiplied  by  appropriate  scalar 
func  tions . 

We  must,  therefore,  study  how  we  can  conveniently  specialize  V.  The 
following  theorem  begins  this  process. 

Theorem  3 — If  U  is  /(-hermi ti an  at  all  z,  then  V  may  be  chosen  to  be 
f(-hermitian  also  at  all  z. 

This  follows  quickly  from  the  fact  that  any  matrix  may  be  separated 
into  two  components: 
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V 

-  V,  +  ;Vj 

(109) 

V, 

■  —  (V  +  K-*V<K) 

2 

(110) 

V, 

=  -  y;(V  -  K"‘VtK) 

(111) 

where  Vj  and  Vj  are  both  /f-hermitian. 

Further,  if  U  is  everywhere  ff-hermitian,  so  is  dV/dz,  as  may  be 
seen  by  differentiating  Eq.  (87),  and  remembering  that  K  is  constant. 
Hence  “j  tV,  Ul  is  A(-hermitian,  so  that 

-;K(V,  +  yVj)!)  +  jKlKVj  +  jV^) 

=  jU»(V/  -  jV2»)K  -  ;(V/  - 
-;KVjl)  +  KVjL  +  jKUVi  -  KUY^ 

=  jU^Vj^K  +  U%»K  -  jV,*U«K  - 
-  ;KUVi  +  KUVj  -  ;KV,U  "  KV^U 

since  U,  Vj ,  and  Vj  are  all  /('hermi tian. 

Hence 

2KVjU  -  2KUV, 


(112) 


(113) 


or  Vj  commutes  with  U. 

But  we  can  add  to  Y  any  matrix  that  commutes  with  U  without  changing 
Eq.  (95)  in  any  essential  way.  Specifically,  we  can  add  “^Vj  and  be  left 
only  with  the  /(-hermitian  Vj . 

We  can,  without  loss  of  generality,  require  V  to  be  everywhere  K- 
hermiti an  i f  1)  is . 

This  constraint  on  Y  implies  a  constraint  on  the  eigenvectors  of  U, 
as  indicated  in  the  following  theorem. 

Theorem  4  —  If  U  is  /(-hermitian  for  all  z,  and  if,  in  Eq.  (95),  V  is 
chosen  to  be  /f-hermitian  for  all  z,  then  the  A(>products  of  the 
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corresponding  eigenvectors  and  generalized  eigenvectors  of  U  are  in¬ 
variant.  In  particular,  it  is  po.ssihle  to  so  choose  them  so  that  they 
are  everywhere  /f-or thogona 1  and  maximally  normalized. 


Consider  the  derivative  of  such  a  product; 


dz 


(u.'Ku^.) 


du  .  ♦ 

— +  u/K - 

dz 


=  -;.i,»(KV  -  V*K)it..  014) 

where  we  have  used  Eq.  (98)  of  Theorem  1.  If  V  is  /(-hermi tian ,  the  right 
hand  side  vanishes  and  we  must  have 


^  Xj.  (115) 

where  is  constant,  independent  of  z. 

On  the  other  hand  since  U  is  /f-hermi tian ,  we  can,  for  any  partic¬ 
ular  value  of  z,  choose  its  eigenvectors  and  generalized  eigenvectors 
so  that 

Jij  =  016) 

as  stated  in  Eq.  (94).  If  the  set  is  chosen  so  that  this  is  true  at  one 
value,  and  if  V  is  K-hermi tian,  we  have  just  shown  that  it  must  be  true 
at  all  values.  Hence  the  set  {u . }  can  be  chosen  to  be  Af-orthogonal  and 
maximally  normalized  at  all  z. 

Conversely,  since  U  is  assumed  everywhere  K-hermi tian ,  we  can  then 
choose  the  set  {u^}  to  be  /(-orthogonal  and  maximally  normalized  everywhere. 
If  we  do  so,  then  the  V  corresponding  to  these  vectors  will  necessarily 
be  K-hermitian  everywhere.  This  follows  by  differentiating  Eq.  (94)  and 
substituting  Eq .  (98).  We  find  that 

Uj»(KV  -  V*K)u,^^.  -  0  (117) 

for  all  i  and  j  and  at  all  z.  .Since  the  set  is  complete,  we  must  have 
V  /(-hermitian  everywhere. 
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This  property  of  /f-orthogonality  and  invariant  maximal  normalization 
permits  us  to  write  V  explicitly  as 

du  . 

V  =  J  5:  cr.— .  (118) 

t  dz 

We  can  do  this  formally  by  postroultiply ing  Eq .  (14)  by  and  adding 

over  all  i.  Sincci  for  a  /f-orthogonal  and  maximally  normalized  set 

5:  o-.u.u^.tK  -  I  (119) 

t 

Equation  (118)  follows. 

Since  V  is  K-hermitian,  we  can,  alternatively,  write 

du  .  ^ 

V  =  K-'VtK  -  -j  2  a.u^.—^K  (120) 

or,  perhaps  better  since  it  explicitly  shows  V  to  be  A(-hermitian , 

We  should  also  note  that  we  are  still  very  far  from  a  complete 
specification  of  either  the  set  .{u.}  or  of  V.  Regarding  V,  we  can  still 
add  to  it  any  matrix  function  of  z  that  is  K-hermitian  and  that  commutes 
with  U.  In  particular  we  can  add  any  suitable  real  function  of  U. 

Regarding  the  set  {u^.},  we  can  make  a  substitution  of  the  form  of 
Eq.  (106),  where  the  are  scalar  functions  of  z,  without  altering 
either  the  orthogonality  or  normalization  properties  providing  only  that 

»  4>i*  .  (122) 

We  have,  still,  an  enormous  amount  of  freedom  available  to  us.  We 
would  like  to  be  able  to  use  this  remaining  freedom  to  assure  that  the 
resultant  V  is  expressible  by  an  equation  similar  to  Eq.  (95) — i.e,,  so 
that  the  resultant  V  has  constant  eigenvalues.  We  have  not  so  far  found 
a  general  way  to  do  this. 
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We  can  push  the  problem  somewhat  further,  however.  If  we  differ¬ 
entiate  Eq.  (121),  we  obtain 


d,  ■  2  '  ^  -< 


u,v,.  K  -  u, - K  I 


If,  now,  it  is  possible  to  define  a  matrix  W  so  that 


dz^ 


=  -2jW- 


dz 


(123) 


(124) 


where  W  is  /(-hermitian,  then  we  can  substitute  this  in  Eq.  (123)  and  use 
Eqs.  (118)  and  (120)  to  give 


dz 


-j[w,  V] 


(125) 


The  problem,  therefore,  is  to  so  specify  the  {u.}  after  a  transfor¬ 
mation  of  the  form  of  Eq.  (106), so  that  W,  defined  by  Eq.  (124)  is  K- 
hermitian.  This  is  a  problem  that  will  require  further  study. 

This,  then  completes  the  general  development  of  the  technique  of 
analysis  in  its  present  state.  Before  considering  a  specific  example, 
we  will  discuss  briefly  the  technique  of  obtaining  solutions  in  those 
cases  in  which  the  sequence  of  Eqs.  (85),  (86),  (90),  etc.,  is  short. 

D.  SOLUTION  OF  THE  SIMPLER  CASES 
1.  R  =  Rq,  Constant 

The  simplest  of  all  systems  of  this  sort  is  when  R  itself  is  con¬ 
stant,  Rj.  This  is  the  uniform  case,  and  has  the  formally  simple  solution 

x(z)  =  e  ^*'’'x(0)  .  (126) 

There  may  be  a  practical  difficulty  in  obtaining  the  exponential  in 
convenient  form,  but,  at  least  in  principle,  it  can  always  be  done. 
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2.  A  =  Aj,  Constant 

The  next  case  is  where  the  A  of  Eq.  (86)  is  constant.  The  solution 
to  Eq.  (86)  is,  then 


It 


;  A, 


(127) 


where  Rq  is  R(0),  Again  there  is  no  difficulty,  at  least  in  principle, 
since  A^  commutes  with  the  exponential,  and  the  derivative  can  be  taken 
without  difficulty. 

To  solve  Eq.  (85)„  then,  we  let 


x(z)  =  e  ^  ®  y(z) 


so  that  Eq.  (85)  becomes 


■je 


Any 


-jAgidy  _  _ 


dz 


= 


•iAfl. 


Roy 


(128) 


or 


dz 


-;(Ro  -  Ao)y 


(129) 


and 


“I* 

y(z)  ■  c  ®  y(0)  >  (130) 

Hence  we  find  that 

x(z)  =  e  ®  e  ®  °  x(0)  .  (131) 

[This  form  is  different  from  that  given  in  Part  II-A  of  Interim  Report  2, 
Eq.  (97).  It  may  he  shown  to  be  equivalent,  liowever,  and  appears  to  be 
a  more  convenient  form.} 
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3.  B  =  Bq,  Constant 

The  next  simplest  case  is  when  the  B  of  Eq.  (90)  is 

As  above,  we  have,  then. 


If,  in  Eq.  (86),  we  let 


R 


and  note  that  D^v  commutes  with  its  exponential  so  that 


dz 


“jB^c 


Then  Eq.  (86)  becomes 


dS 

dz 


.  -j(B  -  AJS  +  ;S(A.  -  BJ 


Hence  we  find  that 


and 


R 


-iBfl* 

e  e 


y(Ao-Bp),|^ 


(A„-Bp), 


If,  now,  in  Eq.  (85)  we  let 


X 


-/B, 


)i 

y 


we  obtain 


constant , 


(132) 


(133) 


(134) 


(135) 


(136) 


(137) 
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dz 


~J^oy 


je 


AflC 


(138) 


This  expression  can  be  solved  by  the  lemma  that,  if  M  and  N  are 
constant  matrices,  then 


e'“Ne“«*  =  I  [  M,  N]zVn! 

it  =  0  " 


(139) 


where 


[*M.  N]  is  the  ^-commutator 


that  is,  where 


[„M.  N]  =  N 


(iM,  IVJ  =  rM,  IV] 


IM,  N]  *  [M[M,  Nl] 


+  N]  =  [M[*M.  NJ] 


(140) 


lienee 


dy 


”  ~jKy  ^  Moy  ^  j  "  ®o)-  Ao]z"/n!)y 


=  "/<Ro  -  Ao)  +  J  f  c, -^y 

n  -  1  n  I 


(141) 


where 


C„  =  [„>(A„  -  Bo).  AoJ  (142) 
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Tlieti 


y 


^  fO 


(143) 


and  X  is  obtained  as 


■jBgi  -j(Ag-B„)T  ->{(R„-Ap).-J.C^2"*'/(n+l  )!> 


(144) 


The  complexity  of  this  solution  is  formidable..  Nevertheless  it 
is  an  exact  solution  to  a  problem  that  admits  a  considerable  degree 
of  non-uniformity. 

We  shall  not  attempt,  at  this  time,  the  exact  solution  of  the 
next  higher  class  of  problems. 


E.  NON- UNIFORM  TRANSMI .SSION  LINE 

We  shall,  now,  consider  a  specific  example:  a  transmission  line 
whose  impedance  is  a  specific  function  of  z. 

This  problem  is  one  of  general  importance,  primarily  because  of 
the  paucity  of  known  soluble  cases  in  spite  of  the  importance  to  micro- 
wave  engineering  of  various  impedance-transforming  or  matching 
transmission-line  sections.  Our  concern  here,  however,  is  to  determine 
the  significance  of  the  various  operations  we  have  performed,  and  to 
gain  insight  into  the  nature  of  the  physical  problems  to  which  this 
analytic  technique  applies. 

In  a  lossless  transmission  line  (single  mode),  described  in  the 
E-I  basis,  the  ||  matrix  is 


(■ 

\/^/^  0  ) 


(145) 


where  [i  is  the  propagation  constant  and  7,  is  the  characteristic 
equation. 

We  shall  assume  /J  to  be  constant  and  allow  7,  to  vary  with  z. 
Since  the  eigenvalues  of  R  are  ±/l,  this  puts  us  in  position  to  apply 
the  techniques  developed  here. 
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It  may  be  seen  that  a  suitable  metric  for  this  system,  independent 
of  2 ,  is 

which  may  be  found  to  lead  to  the  conservation  of  energy. 

We  could  consider  the  eigenvectors  of  Kq.  (145).  However,  we 
might  as  well  consider  a  somewhat  more  general  situation.  It  is  not 
difficult  to  determine  that  the  general  matrix  M(2)  which  is  Af-hermitian 
and  which  has  constant  eigenvalues  is  given  by 


M(z) 


(147) 


where  a  is  a  real  constant  and  /,  g,  h  are  real  functions  of  2  such 
that 


~  gh  -  b  ~  =  c 


(148) 


Here,  b  is  a  real  constant..  (The  constancy  of  the  eigenvalues  requires 
that,  in  the  characteristic  equation ,  the  coefficient  of  each  power  of 
k  shall  be  constant.) 


We  could,  now,  proceed  by  finding  the  eigenvectors  and  going 
through  the  full  process.  In  the  present  case  involving  only  2^2  matrix, 
it  is  easier  to  use  brute  force.  We  seek  an  N  of  the  form  similar  to 
Kq.  (147) 


w 


'/I  +  jF  G 
II  A  -  JF 


(149) 


such  that 


</2 


-j(,N,  M] 


(150) 
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Wc  fitifl  that  we  require 


/'  =  lie  -  Gh 

g'  =  2Fg  -  2Gh 

h'  =  211  f  -  2Fh  .  (151) 

These  equations,,  together  with  the  analogue  of  Kq.  (146) 

-  Gll  =  C  (152) 

where  C  is  a  real  constant,  determine  the  possible  forms  of  F,  G,  and  //. 
They  are  consistent  since,  if  we  substitute  them  in  the  derivative  of 

Eq.  (146) 


2//'  -  gh'  -he'  =  0 

we  do  find  that  the  equation  vanishes  identically,  ' 

We  could  solve  Eq.  (151)  in  terms  of  any  one  of  the  functions. 
However,  to  do  so  would  imply  the  non-vanishing  of  various  terms,  which 
we  would  rather  not  do. 

Let  us  consider  now,  our  original  problem.  To  put  the  R  of 
of  Eq.  (145)  in  the  form  of  M>  we  set 

g  = 

h  =  /i/Z 

/  =  0  (153) 

If  we  substitute  these  into  Eq.  (151)  wc  find  that  we  can  write; 


(154) 
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where 


7/ 

~  2Z 


—  In 
dz 


(155) 


and  tv  is  an  arbitrary  real  function  of  z.  The  condition  of  liq.  (152) 
requires  that  tv  be  of  the  form 


w 


(156) 


where  c  is  an  arbitrary  constant,  restricted  only  to  being  not  greater 
than  the  minimum  of  . 

The  matrix  A  is  constant  only  if  k  is  constant  and  we  set  c  =  , 

By  Eq.  (155)  this  is  the  case  of  an  exponentially  tapered  impedance, 
which  is  too  well  known  to  be  interesting. 

If,  now,  we  consider  the  next  order  of  complexity,  we  find  that 
we  must  satisfy 


ui(ZH)  ~  IV 


IV*  +  2iv^ 


~2wF  +  2x(ZH)  -  v>'  ~  2wx 


(157) 


These  equations  are  consistent.  If  we  subtract  2k  times  the  first 
equation  from  v  times  the  sum  of  the  second  and  third,  we  are  left 
with 


KK '  =  IVtV* 


(158) 


and  this  is  true  by  Eq.  (156)  which  holds. 

The  simplest  solution  to  Eq.  (157)  is  obtained  by  setting  G  “  0. 

Then 


// 


wZ 


F 


t 

w 

=  -  K 

2v) 


(159) 
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For  B  to  be  constant,  F  must  be  constant.  Substituting  from 
Fq,  (156)  for  w' ,  we  must  have 


kk' 

-  +  K 

-  c) 


$ 

c 


or 

KdK 

-  B 

(c*  -  K)(c  -  K^) 


(160) 


or, 


(c  -  c*^) 


In  (c  ~  k)  ~ 


2(c'^  +  c') 


In  (c'^+  K) 


2(c^^-  c') 


In  ~  k)  «  a  -  z 


(161) 


If,  on  the  other  hand,  we  set  =  c' ,  then  the  solution  is 


1/2  ,  1  -  K 

(c'  -  K)  4c'  "  c'  +  /< 


a  ~  z 


(162) 


The  required  variations  of  impedance  are  found  by  integrating  k 
to  gi  ve  In  . 

Other  solutions  are  possible,  which  remain  to  be  worked  out..  These 
solutions  are  sufficient,,  however,  to  demonstrate  that  there  is  a  class 
of  impedance  variations  that  are  soluble,  since  they  give  a  constant 

B  in  Eq .  (90). 

We  have  shown  that  there  does  exist  a  class  of  non-uniform  systems 
that  ore  exactly  soluble.  We  have  shown  some  of  the  properties  of 
these  systems,  and  have  opened  the  way  to  their  more  detailed  analysis. 

There  remain  some  serious  questions  of  a  fundamental  nature..  These 
were  outlined  in  the  introductory  section.  We  have  not  yet  explored 
the  specific  case  of  a  non-uniform  transmission  line  in  sufficient 
depth  to  give  an  adequate  '* feel ”  of  the  significance  of  the  various 
si tuat  ions. 
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Thus  tlicro  ri'iiuiitis  a  cnti  s  i  il(!  I’a  1)  I  o  ainoiinl.  ol  additional  i  ii  vos  t  i  ga  t  i  on 
that  is  noeded.  ilowovcr,  the  approach  docs  seem  tf»  he  important  in  that 
it  does  provide  means  for  e  I  tie  i  da  t  i  iif;  some  of  the  phenomena  that  are 
possible  in  non-uniform  systems. 


I 


I 
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Ill  PROGRAKl  FDR  THE  NEXT  INTERVAL 


Durinfj  tlie  next,  interval  we  will  continue  work  on  the  following 
items; 


(1)  The  multiple  pumping  situation  described  in  Sec.  II-A  of  the 
First  Interim  Report  will  be  considered  further.  We  will 
attempt  to  determine,  at  least  in  our  own  minds,  if  this 
leads  to  a  practical  design  of  a  potentially  useful  device. 

(2)  The  analysis  of  commutator-derived  systems  given  in  Sec.  II 
will  be  continued.  In  particular,  we  shall  look  for  ways 
to  determine  the  gross  properties  of  such  systems  so  that 
the  implications  of  this  analysis  can  be  interpreted  in 
terms  of  physical  behavior  and  device  possibilities. 

(3)  The  First  Technical  Note  on  the  topological  analysis  of 
pairwise  coupled  systems  will  be  issued. 

(4)  Further  study  will  be  made  of  the  group  theoretic  approach, 
if  time  permits,  to  determine  what  use  can  be  made  of  it. 

(5)  Surveillance  of  the  current  literature  will  be  continued. 


j 

i 
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APPENDIX 


THE  ELIMINATION  OF  SELF-GOUPLING  TERMS  IN  R 


52 


APPENDIX 


THE  ELIMINATION  OF  SELF-COUPLING  TERMS  IN  R 


We  shall  here  show  that  it  is  always  possible  to  remove  the  self- 
coupling  terms  from  a  representation  of  the  system  operator.  That  is, 
if  the  system  is  described  by 

=  -;Rx  (A-1) 

dz 

and  R  is  expanded  in  terms  of  a  set  of  constant  A-orthogonal  maximally 
normalized  vectors  {u^} 


R  •  r  .^.cr.u.u.^^.*K  (A.2) 

then  it  is  always  possible  to  eliminate  the  terms  in  ♦K.  We  call 

these  the  '* sel f -coupling"  terms  since  they  produce  a  change  in  the  a. 
component  of  x  that  is  in  the  direction  of  a.. 

To  simplify  the  symbology  somewhat,  let  us  set 

Ej^  -  cr.u.a^,/K  (A-3) 

It  will  be  noted  that  the  E_  are,  then,  idempotent,  and  the  E^yCy/i) 
are  nilpotent.  More  generally: 

E.A*  =  .  (A.4) 

Then  Eq,  (A-2)  can  be  written 

R  =  >;  r..(z)E.  .  (A-5) 

1'he  elimination  from  R  of  the  self-coupling  terms  then  depends  on 
the  following  lemma; 
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Lemma  i:  If  is  a  scalar  and  and  idempotent  dyad,  then 


e 

*  _  j 

-  (1  -  e 

This  is  proved  by 

the  series  expansion  of  the 

idempotency  of  j 

e  = 

I 

l*k^kk 

—  V2 
2!  ** 

+  ...  T(- 

n 

= 

I 

/ 

+ 

1  -hn,  - 

2! 

c- 

u-jr  — 

n  n ! 

>“k 

“ 

I 

“  E**  +  e 

Ekk 

s 

I 

-  (1  - 

(A-6) 


n! 


F" 

I 


As  an  aside  which  we  sliall  not  use  liere,  but  for  possible  future 
reference,  we  may  note  the  following  lemma: 

Lemma  2:  If  't  is  a  scalar  and  a  nilpotent  matrix,  then 


e 


aE  . 

'  *7 


I  -  >«E, . 


(A-7) 


Again,  we  prove  tlie  lemma  by  expanding  the  exponential  in  a  power 
series,  this  time  obtaining  the  results  directly. 

Returning  now  to  our  immediate  purpose,  consider  a  change  of 
variable  given  by 


X 


(A-S) 


where  is  a  function  of  z  to  be  chosen  as  desired.  We  may  note  that 
tills  change  affects  only  the  component  of  x.  If  ^  i,  '  s  real  it  amounts 
to  making  the  reference  f»hase  angle,  agninst  whicli  the  phase  of  this 
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component  is  measured,  a  function  of  position.  If  tlie  is  complex, 
then  there  is  in  addition,  a  z>dependent  renormalization  of  the 
component  of  x. 

We  may  substitute  this  into  Eq.  (A-1)  with  Eq.  (A-5)  for  R  and 
obtaining  the  derivative  of  the  exponential  from  Eq.  (A-6) 


dct 

dz 


k  - 
e 


„  +  ~’°'k^kk 


dy 

dz 


“7  2 


r  .  E.  .{I  -  (1  “ 

t  j  t  j 


- 1  a . 

C  *)E 


kk- 


1 

i 


I  k 


E 


ikY 


r . .E .  .y 


(A-9) 


If,  now,  we  set 


(Xj  =  J  r^^[z)dz  (A-10) 

then  the  first  term  on  the  left  cancels  the  term  in  the  first  sum  on  the 
right  with  t  ■  A.  There  remains  the  following  equation: 


^ijaktkkjy 


■  K- 

-}  ^  r.jC  ‘Ej* 
I* 


dz 


L  . 

tfk  ' 


(A-in 
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I’remu  1 1  i  p  1  y  i  ng  by 


dz 


If,  now,  we  set 


then  we  have 


by  comparison  to  Kq 


,  cxpaniJi'il  by  I'lq.  (A-6),  we  obtain 


-j(l  -  (1  - 


-j(I  -  (1  - 


.  •C- 

-J  r  r.  .e  *E.  .y 


j  y.  r  .  E.  .y 

jfk 


(A-12) 


-  0  if  i  =  ;  =  t 

=  r.^. e  if  i  /  k,  j  =  k 

=  ^  if  i  =  fe,  ;•  /  k 

•=  r.^.  if  i  ^  k,  j  i  k  (A-13) 

I77  =  ■>-'’>  "  "■'"'y  <A.14) 

.  (A-5). 
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We  have,  in  fact,  eliminated  from  R  the  term  in  —  i ,  e  .  ,  the 
self-coupling  feth  component  of  R.  The  other  terms  are  multiplied  by 
scalar  factors  if  they  couple  to  or  from  (i  “  k,  or  j  “  k) ,  and  are 
unaffected  otherwise.  In  particular,  the  other  se 1 f -coupl ing  terms, 
i.e,,  those  in  f  k) ,  are  unaffected. 

We  can,  therefore,  proceed  in  the  same  manner,,  with  the  elimination 
of  the  term  in  without  reintroducing  a  term  in  E^j^. 

Since  the  dyads  E^  j  commute,  the  law  of  combination  of  exponentials 
applies.  We  can,  therefore,  accomplish  the  simultaneous  elimination  of 
the' sel f-coupled  terms  by  the  substitution 

X  =  exp  -  J  a^E^y 

where 

We  have,  therefore,  a  procedure  for  the  complete  elimination  from  R 
of  all  self-coupling  terms. 
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1.  K.  Blotekjer,  "Transverse  Electron  Beam  Noise  Described  by 
Filamentary  Beam  Parameters,”  J.  Appl.  Phys.  33,  2409-2414 
(August  1962). 

Summary  —  The  effect  of  transverse  emission  velocities  and 
positional  variation  (transverse  shot  noise)  on  an  equivalent 
filamentary  beam  is  studied.  The  frequency  dependence  of  the 
noise  in  the  various  beam  modes  is  obtained.  It  is  shown  that 
there  is  an  effect  which  is  similar  to  the  space-charge  smoothing 
of  longitudinal  shot  noise  which  acts  to  reduce  the  effect  of 
positional  variation.  This  effect,  however,  does  not  seem 
sufficient  to  account  for  the  low  noise  observed  in,  for  example, 
the  Adler  tube.  Hence  the  author  concludes  that  "successful 
operation  of  the  eye lotron -wave  parametric  amplifier  is  largely 
due  to  an  effect  which  is  not  yet  understood.  ” 

Comment — A  useful  analysis  of  the  generation  of  the  noise 
that  is  the  input  to  devices  using  various  mode-coupling..  Of 
interest  also  is  the  clear  result  that  existing  theories  are  not 
wholly  adequate. 

2.  D.  C.  Forster,  “Cooling  of  the  Slow  Space-Charge  Wave  with 
Application  to  the  TWT.  ”  I/?£  Trans.  PGED-9,  449-453  (November  1962). 

.Summary — An  analysis  of  the  beam-cooling  method  originally 
proposed  by  Sturrock.  In  this  type  of  device,  a  fast  wave  is 
first  cooled  by  a  Kompfner-null  coupler.  A  slow  wave  is  then 
coupled  to  the  fast  wave  parametrically.  Because  of  the  parity 
inversion  under  parametric  coupling,,  the  resultant  coupling  is 
passive,  or,  in  our  terminology,  /3-coupling,  Hence  the  noise  on 
the  slow  wave  is  transferred  to  the  fast  wave..  The  resultant 
cooled  slow  wave  can  then  be  used  in  a  conventional  TWT  structure. 
The  analysis  given  here  is  based  on  linear  coupled-mode  theory, 
and  involves  a  detailed  study  of  the  effects  of  the  various  coupling 
terms  that  might  be  expected.  The  results  of  a  computer  calibration 
are  presented,  and  indicate  a  possible  noise  temperature  of  about 

lOO^K. 
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Comment — As  far  as  wc  know,  the  feasibility  of  beam  cooling  by 
this  method  has  not  yet  been  demonstrated  experimentally.  This 
paper  demonstrates  to  a  reasonable  degree  that  the  difficulty,  if 
there  is  one,  is  due  to  the  neglect  of  higher-order  linear  terms. 
What  is  perhaps  needed,  instead,  is  the  investigation  of  the  effect 
in  the  much  more  complicated  mode  structure  found  recently  by 
Tore  Wessel-Berg  by  studying  thick  beams.  In  other  words,  this 
paper  is  interesting,  but  may  be  considering  the  wrong  problem. 

3,  4.  D.  L.  Uobroff,  II.  A.  Haus,  and  J.  W.  Kliiver,  "On  the  Small  Signal 
Power  Theorems  of  Electron  Beams,"  J,  Appl.  Phys,  33,  2832  (1962), 
and  J.  W,  Kliiver  “Potential  Form  of  the  Small  Signal  Power  Theorem," 
J.  Appl.  Phys.  33,  2943  (1962). 

Summary — The  small  signal  power  theorem  for  electron  beams 
has  been  a  subject  of  controversy  for  about  three  years..  Certain 
difficulties  in  the  previous  works  by  the  authors*'*  on  this  subject 
have  been  discussed  by  E.  L.  Chu.**®  The  authors  have  subsequently 
answered  in  rebuttal.* 

The  present  papers  do  not  appear  to  clarify  the  situation. 
Little  material  not  previously  published  is  included,  although  the 
mathematical  form  of  the  power  theorem  has  been  favorably  modified.. 
Instructive  examples  are  presented. 

The  following  is  a  bibliography  of  some  of  the  pertinent 
articles  on  this  subject: 

1.  H..  A.  Haus  and  D.  L.  Bobroff,  “Small  Signal  Power  Theorem 
for  Electron  Beams,"  J.  Appl.  Phys.  28,  694  (1957). 

2.  J.  W.  Kliiver,  "Small  Signal  Power  Conservation  Theorem  for 
Irrotational  Electron  Beams,  "7.  Appl.  Phys.  29,  618  (19.58). 

3.  E.  L..  Chu,  “Two  Alternative  Definitions  of  Small  Signal  RF 
Power  of  Electron  Beams,"  J.  Appl.  Phys.  30,  1617  (19.59). 

4.  E.  L.  Chu,  "Comments  on  Kliiver' s  Paper  Entitled  Small  Power 
Conservation  Theorem  for  Irrotational  Electron  Beams," 

J.  Appl..  Phys.  30,  1618  (1959). 

5.  E.  L.  Chu,  “The  Lagrangian  and  thy;  Energy-Momentum  Tensors 
in  the  Perturbation  Theory  of  Classical  Electrodynamics," 

Annals  of  Phys,  9,  76  (I960), 


59 


6.  i;.  I..  Chu,  “On  tlie  (Concept  of  Kictitious  Surface  Charges  of 
an  KJectron  Ueatn,”  J.  Appl.  Phys!  31,  381  (1960). 

7.  D.  L.  Bobroff,  II.  A.  Ilau.s,  and  .1.  W.  Kluver,"0n  E.  L.  Chu’s 
Definition  of  Sma 1 1 -Si gna 1  DF  Power  of  Electron  Beams," 

J.  Appl.  Phys.  .32,  749  (1961  ). 


60 


STANFORD 

RESEARCH 

INSTITUTE 


MENLO  PARK 
CALIFORNIA 


Regional  Offices  and  Laboratories 


Southern  California  Laboratories 

820  Mission  Street 
South  Pasadena,  California 

Washington  Office 

808  17th  Street,  N.W. 

Washington  5,  D.C. 

New  York  Office 

270  Park  Avenue,  Room  1770 
New  York  17,  New  York 

Detroit  Office 

The  Stevens  Building 
1025  East  Maple  Road 
Birmingham,  Michigan 

European  Office 
Pelikanstrasse  37 
Zurich  1,  Switzerland 

Japan  Office 

911  lino  Building 

22,  2-chome,  Uchisaiwai-cho,  Chiyoda-ku 
Tokyo,  Japan 


Representatives 


Honolulu,  Hawaii 

Finance  Factors  Building 
195  South  King  Street 
Honolulu,  Hawaii 

London,  England 

19  Upper  Brook  Street 
London,  W.  1,  England 


Milan,  Italy 

Via  Macedonio  Mellon!  40 
Milano,  Italy 

London,  Ontario,  Canada 
P.O.  Box  782 
London,  Ontario,  Canada 


